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ABSTRACT 

Rezzolla et al. [ApJ 531 (2000), L139; Phys. Rev. D 64 (2001), 104013; Phys. Rev. 
D 64 (2001), 104014] draw attention to the second order secular drift associated with 
r-modes and claimed that it should lead to magnetic field enhancement and sup¬ 
pression of r-mode instability in magnetized neutron stars. We critically revise these 
results. We present a particular second order r-mode solution with vanishing secular 
drift, thus refuting a widely believed statement that secular drift is an unavoidable 
feature of r-modes. This non-drifting solution is not affected by magnetic field B, if 
B <C Bcrit ~ 10 17 (^/600Hz) G (u is a spin frequency) and does not lead to sec¬ 
ular evolution of magnetic field. For general second order r-mode solution the drift 
does not necessarily vanish, but the solution can be presented as a superposition of 
two solutions: one describes evolution of differential rotation in nonoscillating star 
(which describes secular drift; for nonmagnetized star it is arbitrary stationary ro¬ 
tation stratified on cylinders; for magnetized star differential rotation evolves on the 
Alfven timescale and may lead to magnetic energy enhancement), and another one 
is non-drifting r-mode solution mentioned above. This representation allows us to 
conclude that enhancement of magnetic field energy is limited by initial energy of dif¬ 
ferential rotation, which is much less (for a factor oc a 2 , where a is mode amplitude) 
than the total energy of r-mode. Hence, magnetic field enhancement by drift cannot 
suppress r-mode instability. Results can be generalized for any oscillation mode in any 
medium, if this mode has non-drifting solution for B = 0. 
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1 INTRODUCTION 

Observations of millisecond pulsars demonstrate that some 
neutron stars (NSs) rotate very rapidly (the fastest known 
pulsar P SR J17 48-244 6ad has spin frequency v — fl/27r ~ 
716 Hz; Hessel s et al.l 2006|) . However, as was s hown by 
lAndersson! ( 19981 ) and iFriedman fc Morsinkl d 19981 1. rapidly 
rotating NSs are subject to gravitation-driven instabil¬ 
ity associated with enhancement of r-modes (toroidal 
oscillation mode of rotating star controlled by Corio¬ 
lis force). It is a particular case of the Chandras ekhar- 
Friedman-Schutz (CFS) i nstability [Chandrasekhar Il97d : 
IFriedman fe Schut d 1 1978a .bl. R-mode instability probably 
plays an important role in NS physics. Indeed, it is ar¬ 
gued to be a m e chani s m which limits pulsars’ spin frequen ¬ 
cies |Bildgten[_Jl99^jAj2dnrssoj2 ^Kokkotas fe Stergioulasl 
Il999l : iGusakov, Chugunov fe Kantoil 2014bl a). Rapidly ro¬ 


tating hot NSs, observed in low mass X-ray bi¬ 
naries (LMXBs), can be affected by r-mode insta- 
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bility even more profoundly: some of these stars 
should be unstable wit hin the standard model of r- 
rnode oscillations (see w. i Andersson_&Haskel| 1201 ll : 


iHaskelL Deeenaar fe Holl2012l ; Gusakov et al. 2014bl . for ex¬ 
ample) and their existence provides one with a good oppor- 
tunity to test microphysical models of NS’s core (see, e .g., 
IHo et al.ll201^ Haskell_et_ah 201^ Gugakov etalj 201 Tbl fal) . 


Furthermore. Ichueunov. Gusakov fe Kantoil 1 20141 1 suggest 

that r-mode instability can support a new class of rapidly 
rotating NSs, which do not accrete matter from com¬ 
panion even transiently (as NSs in LMXBs), but still 
keep high temperature via r-mode instability. Some of 
these NSs might have been already observed, but er- 
roneously classified as quiescent LMXB candidates (see 
IChueunov, Gusakov fe K antoi||2014|fo r deta iled discussion). 
Finally, as argued by iRezzoIla et all (|2000l l. r-mode insta¬ 
bility can be also important for generation of NS magnetic 
fields. The main aim of the paper is to discuss the role, which 
magnetic field plays in nonlinear evolution of r-modes and 
role of r-modes in magnetic field evolution in NSs. 

Most of the papers discussing r-modes in NSs either 
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2 Audrey I. Chugunov 


neglect magnetic field or simply mention that it can be 
important. Indeed, dipolar magnetic field of rapidly rotat- 


ing NSs (u > 200 Hz) is not very large ~ 10 s G and 
seems to be clearly negligible (see Rezzolla, Lamb & ShaDird 

20001: Asai, Lee & Yoshida 

2015 and Sec. (3l). Recent 

papers (Lee 2005 

: Abbassi. Rieutord & Rezania 

20 lj; 

Chirenti & Skakala 

20131; Asai, Lee & Yoshida 2015 

, for in- 


stance) confirm by detailed numerical calculations that only 
very high magnetic fi eld can affect r-mod e s stro ngly [for ex¬ 
ample (see fig. 7 by IChirenti fe Skakalal l|2013h l. magnetic 
field B > 3 x 10 15 G modifies r-mode frequency for a less 
than for a one percent for v ~ 220 Hz], H owever, just af¬ 
ter r-m ode instability was discovered, iRezzolla et al.l d2000l . 
l2001alBh pointed that r-modes are able to increase seeding 
magnetic field through secular drift of fluid elements, which 
takes place within second order (in oscillation amplitude) 
perturbation theory. They argue that this mechanism can 
be responsible for generation of strong magnetic fields in 
NSs, but at the same time it finally suppresses r-mode in¬ 
stability. A strong support for these results was provided 
by subsequent papers, which claim differential rotation, i.e. 
secular drift of fluid elements, to be a necessary feature 
of r-mo des (in nonmagnetized NS) on the base of analyt¬ 
ical (Islfcool K fe Tomdl2005ll and numerica l calculations 
(see IStergioulas fe Font 200 ll : iLindblom et alJ l200ll . for ex¬ 
ample). 

In this paper we reanalyze these results using analyt¬ 
ical derivations within second order (in oscillation ampli¬ 
tude) perturbation theory. In Sec. [2] we discuss nonmagne- 
tized ( B = 0) NSs, f ollowi ng the approach suggested bvlsll 
ll2004) ; ISa fe Tome! (120051 1. As they have shown, the sec¬ 
ond order r-mode solution is not unique, but determined 
up to any differential rotation, which can exist in nonoscil¬ 
lating star (i.e. arbitrary stationary differential rotation, 
stratified on the cylinders). However, an important point 
of these papers, namely that “differential rotation, produc¬ 
ing large scale drifts of fluid elements along stellar lati¬ 
tudes, is an unavoidable fea ture of r-modes in the nonlin¬ 
ear theory” (see abstract in |S4 fe Tomell2005lj is not true: 
we present a paramete r set f or analytical secon d order so¬ 
lution obtained byl! d2004h : ISa fe Tomel d2005h with van¬ 
ishing secular drift (below we refer to it as to “non-drifting 
solution” )Q We argue that existence of second order solu¬ 
tion with vanishing drift is not coincidental feature, spe¬ 
cific for r-modes only, but rather is a general property of 
oscillation modes. In Sec. [3] we reconsider the effect of mag¬ 
netic field on r-modes in second order in amplitude by ac¬ 
curate treatment of mag netic stresses w i thin ideal ma gneto- 
hydrodynamics fMHDL IRezzolla et alJ (l2000l . l2001al bh ar- 


1 Isil 120041 ) and lSa &; Toim-j f 20051 ) base their conclusion on the 
oscillation averaged second order perturbations of Eulerian ve¬ 
locity. Indeed, these perturbations are nonvanishing for all possi¬ 
ble second order solutions. However, the drift of fluid elements, 
i.e. secular increase of the Lagrangian displacement, is not deter¬ 
mined by the second order Eulerian velocity only, but has correc¬ 
tion associated with Stokes drift induced by the first order motion 
of the fluid elements (see [Longuet-Higg ins 1953; Rc zzolla et al.l 
l2001al or Eqs. (28) and (T2) in Isal1 20041 and lSa &; Tomdl2005l , re- 
spectively). This correction allows to have second order solution 
with vanishing secular drift, but nonvanishing oscillation averaged 
Eulerian velocity (see Sec. [2] for details). 


gued that r-mode velocity profile is not affected by mag¬ 
netic field up to B 10 16 G. Assuming secular drift to be 
unaffected, they conclude that r-modes enhance magnetic 
field due to differential rotation, associated with drift □ We 
demonstrate that non-drifting solution is indeed unaffected 
by magnetic field up to B <C B cr it « 10 17 (^/600 Hz) G, 
but it does not lead to secular evolution of magnetic field. 
However, in a general second order solution the second order 
velocity profile stays unaffected by magnetic field only for a 
timescale much less then the Alfven timescale (while back- 
reaction of magnetic field can be neglected). On the contrary, 
at the Alfven timescale the drift contribution to the second 
order velocity profile (corresponding to arbitrary station¬ 
ary differential rotation) is modified by magnetic field [back- 
reaction is crucial and cannot be neglected for the same rea¬ 
sons as why stationary differential rotation is forbidden in 
magnetized stars (Ferraro’s law of isorotation)]. To describe 
back-reaction properly, we demonstrate that general second 
order r-mode solution in a magnetized NS with B <C Bcrit 
can be presented as a superposition of two independent solu¬ 
tions of MHD equations: (a) solution which describes evolu¬ 
tion of differential rotation in nonoscillating magnetized NS 
(drift solution, which describes evolution of secular drift), 
and (b) non-drifting r-mode solution. These solutions are 
decoupled (at second order in oscillation amplitude) as far 
as Eulerian perturbations of all quantities (except magnetic 
field) in drift solution are second order in amplitude of non¬ 
drifting r-mode solution. Evolution of the drift solution is 
governed by the same equations as in non-oscillatin g mag¬ 
netiz ed NS, thus it evolves at Alfven timescale fsee ISpruitl 
119991. for example). Secular evolution of the magnetic field 
is solely determined by the drift, thus initial energy of dif¬ 
ferential rotation provides an upper bound for the increase 
of magnetic field energy. Furthermore, energy of the non¬ 
drifting mode is conserved in the leading order. As a result, 
r-modes cannot convert their energy to magnetic field. Hence 
we conclude, that magnetic energy enhancement by secular 
drift cannot suppress r-mode instability^ 


2 R-MODES IN NONMAGNETIZED NS 

In this section we consider r-modes oscillations in 
slowly rotating Newtonian NSs with barotropic equa¬ 
tion of state (i.e. pressure depends only on density). 


2 As it is stressed by I Alfven Sz Felthammar] dl963f) (Sec. 3.14.5), 
the twist of magnetic field caused by inhomogeneous rotation, 
propagates along magnetic lines of force with Alfven velocity ua, 
thus if rotation inhomogeneity is too small, the force lines will 
have enough time to straighten up. As a result, to produce sig¬ 
nificant field by differential rotation, relative velocity of different 
parts of the star should be at least comparable with Alfven ve¬ 
locity. Howev er, this point was not discussed by IRezzolla et al.l 
(2000, 2001a b). On the contrary, they suppose that arbitrary 
small secular dr ift can enhance magnetic field. 

3 As argued bv lMendefll (l200llh [Kinney fc Mendefil (1200315 , mag¬ 
netic field can affect r-mode damping at the crust-core boundary, 
but it requires rather strong radial field (> 10 11 G), which is 
much larger than typical dipolar magnetic fields of rapidly rotat¬ 
ing NSs (« 10 8 G). It is also unlikely that strong radial magnetic 
fields can be generated by drift. 
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R-modes in magnetized neutron stars 3 


In case of general relativity o scillation equations be¬ 
come much more complicated (se e Ruoff fe Kokkotasll2002l ; 


by (also to leading order in fi) 


iLockitch . Friedman fe Anderssonl 120031. for examp le). Fur 


Lockitch. Andersson fe Friedmanl (j200Cli 'l demon 


thermore, 

strate fully relativistic equations for r-modes in NSs with 
barotropic and nonbarotropic equations of state to be qual¬ 
itatively different. For barotropic equation of state a rel¬ 
ativistic analogue for Newtonian r-modes was found by 
iLockitch et al.l (120031 ') and they conclude that “unstable r- 
modes remain essentially unaltered when the problem is 
studied in full general relativity” Q It allows us to take ad¬ 
vantage of Newtonian r-modes for the sake of simplicity. 

Second order (in oscilla tion amplitude) r-mode solution 
was obtained by [si| (120041 '). This solution is applicable if 
oscillation amplitude a (O/Ok) 2 (here Ok ~ \JGM/R 3 
is the mass shedding freq uency ). Several important features 
of r-mode solution by[Saj (12004! ) can be generalized for rather 
general oscillation mode (see below). Thus, below we do not 
restrict our consideration to r-modes, but we apply Sa’s r- 
mode solution as an example. 

Let’s look at the equations describing oscillations up 
to second order in oscillation amplitude a. The first order 
equations are: 

dtS w Vi+8 m v k V k Vi+v k V k S w Vi = -X7i6 w U, (1) 

dtS W p + v i V i S m p + V i {p5 m v i ) = 0, (2) 

A <5 (1) 4> = 4ttG<5 (1) p, (3) 

where v, p, <E>, p are, respectively, velocity, density, grav¬ 
itational potential, and pressure in the unperturbed NS. 
ijWj _ 0(a l ) represents the * th order Eulerian perturba¬ 
tion of quantity /, and S^U = S^p/p + Finally, G is 

gravitational constant. In le ading order in Q , linear r-mode 
solut ion is well known (see IProvost. Berthomieu fe Roccal 
Il98ll . for example) and velocity perturbation can be written 
in forrrQ. 


= i 


f ~lRr 


\/l(l + 1) 


(J)‘vx (rV^)e“ 


(4) 


where Yim is the spherical harmonic with the multipolarity l 
equal to m, Yj m = Yu = yj(21 + l)\/[2 l + 2 it (/!) 2 ] sin l (6)e ll<l> . 
Here uj is the oscillation frequency in the inertial frame, given 


4 IChirenti & Skakalal (l2013ll stated that their relativistic Cowl¬ 
ing calculations give r-mode gravitational radiat ion growth time 
30% l arger than in previous Newtonian results (iLindblom et all 
119981 ; lAndersson et al.lll999h . This statement seems to be par- 
tia lly biased by different radi us and mass of NS taken as fiducial 
by IChire nti fc Skakalal (120131) (R = 14.08 km and M = 1.4M©) 
and by lAndersson et al.l Jl999l) (R = 12.47 kmandA£_=J L: hM^). 
Being rescaled to the fiducial parameter s of | Chirenti & ; Skakala 


2013) by applying scaling relations fro m I An dersso n &; Ko kkotai 
2001 ) [r gr oc 1 /(MR 4 )], results by ILindblom et al.l \ 19981) ; 


lAndersson et ak | Jl999| ) come t o a mu ch better agreement with 
results by IChirenti Sz Skakalal J2013lb Correspondent growth 
time is r gw ~ 13.65(lkHz/z/) 5 - 83 , r gw 12.47(lkHz/z/) 6 , 
13.73(lkHz/z/) 593 , for Refs . [Chirenti fc_S ka.ka.lal i2013h . 


ILindblom et al.l i 19981) , and lAndersson et alj (1 19991) respectively. 


5 We use normalization of the mode amplitude from lOwen et all 
1998), which is usually applied for r-mode instability analysis. I Sal 


2004 h lSa & Toin3 i 20051) apply slightly different normalization, 


<j = _ ( i -i)a + 2) n 

( + i 


(5) 


( 6 ) 


namely their amplitude as a — ol/Y 2. 


The second order (in ck) equations are: 

d t S (2) Vi + 6 {2) v k V k Vi + v k V k S {2) Vi + 5 (1) v k V kl 5 (1) v, 

p V p J 

d t s (2) p + t^Vi^p + ViO^V) 

+ Vi(<5 (1) p5 (1) v i ) = 0, (7) 

A5 (2) $ = 4t rG<5 (2) p. (8) 

Let us stress two features of these equations: 

• For a given first order solution, general solution of sec¬ 
ond order equations is a sum of partial solution of inhomo¬ 
geneous equations [Eqs. J6][8]), which contains terms induced 
by the first order solution] and a general solution of homoge¬ 
neous part of Eqs. (|6lf8|l [without terms induced by the first 
order solution]. 

• Homogeneous part of equations m have the same 
form as linearized equations [Eqs. (IH3I) ]. Thus, evolution of 
all quantities in arbitrary solution of homogeneous equations 
is the same as for (linear) perturbations of initial (nonoscil¬ 
lating) state. 

Thus, with accuracy up to the second order in a, one can se¬ 
lect an arbitrary partial solution of the oscillation equations 
and present a general solution as composition of selected 
oscillating solution and solution, which describes arbitrary 
initial 0(a 2 ) perturbations, evolving in a same way as in 
nonoscillating NS. In fact, it is obvious even without explicit 
equations, since second order perturbations can neither af¬ 
fect oscillation solution nor be affected by it (within a 2 accu¬ 
racy). Let’s restrict subsequent consideration to nonoscillat¬ 
ing solutions of homogeneous Eqs. (I6l8l) . In case of r-modes 
it is stationary differential rotation, stratified on cylinders. 
The freedom in a homogeneous part of the solution can 
be also understood as excitation of r-mode in a differen¬ 
tially rotating star. As far as this background differential 
rotation is weak (second order in oscillation amplitude) it 
does not affect r- mode. Thi s resu lt is in agreement with 
iKarino, Yoshida fe Eriguchil (l200lh . who discuss r-mode os¬ 
cillation of differentially rotating polytropic stars and obtain 
r-mode solution in a case, when degree o f d iffere ntial rota¬ 
tion is not too high. Recently. IC hirenti . Skakala fc Yoshidal 
( j2013l ) ; IChirenti, Skakala fe Yoshida ( 2014 ) confirm that r- 
modes are rather insensitive to small differential rotation 
within general relativity Cowling approximation. Note, how¬ 
ever, that in case of very stron g di fferen tial rotation, r-modes 
can be suppressed (see IKarino et al.l ()200lh . for more de¬ 
tails). 

Partial solution, mentioned above, can have drift, i.e. 
the Lagrangian displacements £ can increase with time pro¬ 
ducing large scale secular motion of fluid elements. Cor- 
responding velocity ( so-called mass transfer velocity; see 
ILonguet-Higgina 11953I . for example) = dt£ is second 

order in oscillation amplitude. Here overline means averag¬ 
ing over oscillations. Generally, mass transfer velocity is not 
equal to oscillation averaged Eulerian velocity, but has a 
correction associated with Stokes drift velocity induced 
by first order pure oscillating solution v^ = (^WV)^Wd 
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4 Audrey I. Chugunov 


(see lLonguet-Higginsil 1953 . for details). As a result, Eulerian 
velocity perturbations can be written in form: 


8^1)' 

= d t ^ 1)i + v j V J ^ 1)i -^ 1)j V j v i 



= (ft + na*)£ (1)i , 

(9) 

8^v l 

= dt^ (2)i + v j Vj^ 2)i 



- ^ 1)j V j S w v i - ^ (1)j £ Wk V j V k v i 



= (d t + nd<t,K (2)i 

(10) 


Here the last equality corresponds to solid state rotation at 
frequency fl for unperturbed state = fl<5^). By averaging 
of these equations one can obtain relations for mass transfer 
velocity. 

For r-modes mass transfer velocit y can be easily 

derived from Eqs. (30-31) in|sij (120041 1. which gives the sec¬ 
ond order Lagrangian displacement. Only (j> component of 
is nonvanishing and it can be written in form0 


(u (mt V = Q 2 n 


A 


r sin# 

R 


21-2 oo 


+ Y, An 


r sin# 

R 


where An are arbitrary constants, and 

(2Z + 1)! (21 - 1 )l 


A = 


221+4. . 


IP 


( 11 ) 

( 12 ) 


For A 21-2 = —A and An = 0 for all other N we obtain 
a second order r-mode solution with vanishing mass trans¬ 
fer velocity^ For this solution Lagrangian displacements do 
not increase with time, thus secular drift of fluid elements 
vanishes. However, oscillation averaged Eulerian velocities 
are n ot zero for this solution (see Eqs. (18) and (21) in Isil 
120041 1. Such a solution corresponds to a very specific set of 
initial conditions and it is not surprisi ng, that it was not 
obtained in numerical calcu lations by IStergioulas &: Fontl 
l|200ll) : lLindblom et al.l (l200ll l. Nonetheless, it simplifies dis¬ 
cussion of r-modes in magnetized NSs (see Sec. [3]). 

Presence of a second order solution with vanishing drift 
is not an accidental specific feature of r-modes. This prop¬ 
erty can be generalized for any oscillation mode, if mass 
transfer velocity for initial partial solution of second 

order Eulerian equations corresponds to a solution of linear 
perturbation equations (or, equally, homogeneous second or¬ 
der equations) 0 I 11 this case one can construct a partial so¬ 
lution with vanishing drift by substituting 8v —» Sv — S^v, 
Sp -+ Sp- S w p, <5<f> Here 8 W p and <5 (d) <f> 

are perturbations of density and gravitational potential in a 


6 In Eqs. (30-31) by ISai l‘2004l i summation over N is omitted, but 
it is required for general solution. 

7 Note, if Ao A 0, it corresponds to uniform rotation, i.e. pertur¬ 
bation of the rotation frequency of the star. Such drift will not 
perturb magnetic field energy. All results obtained for nondrift¬ 
ing solution can be also applied for solutions, which differ from 
nondrifting solution by second order uniform rotation. 

8 It is likely to be rather general case. For example, for ocean 
waves (incompressible fluid in external gravity field) th i s prop erty 
was demonstrated by Longuet-Higgins (1953); Moore (1970). In 
particular, for r-modes this generalization may allow to go beyond 
restriction of second order solution by [si <120041) (n < fl K and 
a ^ (Q/Qk) 2 )- However, we leave formal discussion, whether it 
is a general property or not, beyond the scope of the paper. 


solution of linear perturbation equations, which corresponds 
to velocity perturbation 

Concluding, in this section we demonstrate, that gen¬ 
eral second order r-mode solution can be presented as a 
sum of oscillating solution with vanishing drift and a so¬ 
lution, which describes differential rotation in nonoscillat¬ 
ing NS. As far as we consider nonmagnetized star, this dif¬ 
ferential rotation can be described by arbitrary stationary 
profile, stratified on cylinders. If we restrict our consider¬ 
ation to r-modes excited by the gravitational radiation in 
initially uniformly rotated NSs, this p rofile is well defined 
(lFriedman. Lindblom fe Lockitchl I2015T ). However, station¬ 
ary differential rotation is forbidden for magnetized star 
(Ferraro isorotation law), thus for magnetized NS a general 
r-mode solution cannot be presented as sum of oscillations 
and free stationary differential rotation. Magnetic field mod¬ 
ifies general r-mode solution and we discuss this modification 
in the next section. 


3 R-MODES IN MAGNETIZED NS 

Oscillations in nonmagnetized NS can have nonvanishing 
drift, leading to large (and increasing) Lagrangian displace¬ 
ments, but Eulerian perturbation of all quantities stays 
small. This is not the case for a magnetized NS. Magnetic 
field is frozen into fluid elements, so that large Lagrangian 
displacements generally lead to large Eulerian perturbations 
of magnetic field. As a result, one generally cannot expand 
Eulerian perturbation of magnetic field in powers on oscil¬ 
lation amplitude and apply perturbed Eulerian approach to 
describe r-modes with nonvanishing drift. Hence, a stricter 
formalism is required to describe oscillations with nonvan¬ 
ishing secular drift in magnetized NS 0 Here we describe such 
formalism based on the idea to present a general second or¬ 
der r-mode solution as a superposition of two solutions: (a) 
solution which describes evolution of the differential rotation 
in nonoscillating NS within exact MHD equations (later re¬ 
ferred to as ’drift solution’, because it will describe drift), 
and (b) non-drifting oscillating solution for B = 0, described 
in the previous section. To achieve that, we apply perturba¬ 
tions in two steps. First, we apply initial perturbation, which 
is second order in a, and evolve it using exact MHD equa¬ 
tions. This solution describes evolution of drift motion [i.e. 
it is the drift solution mentioned above]. Second, we per¬ 
turb drift solution, write down perturbation equations up 
to second order in a, and show that non-drifting oscillat¬ 
ing solution of amplitude a is a (partial) solution of these 
equations. Finally, we demonstrate that r-mode with arbi¬ 
trary second order differential rotation at initial moment 
can be presented in form of this superposition, thus the su¬ 
perposition can be applied to describe a general solution in 
magnetized case. It is worth to note that this procedure is 
suitable for any oscillation mode in any media, if the mode 
have non-drifting solution for B = 0. 

Let (p,p, v, B ) be a stationary configuration described 


9 A pure Alfven waves are linear at arbitrary amplitude and one 
can use Eulerian perturbations of magnetic field to describe these 
waves even for large amplitudes. Still, it is not the case of r-modes. 
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R-modes in magnetized neutron stars 5 


by solution of ideal MHD equations^ 


dtv + (vV)v 

= -^-V$+-l[Bx[VxB]], 

p Anp 1 

(13) 

dtp + V(pw) 

= o, 

(14) 

A4> 

= 4nGp, 

(15) 

dB 

= [V x [v x B]]. 

(16) 

dt 


Let c(a, t) be a fluid element trajectory, corresponding to 
this stationary solution^ Solution of the induction equatio n 
can be written in the form (see, e.g.. Iftezzolla et aDl2001al l 

BRc(a,t),t) _ B\a,t = 0) dc J (a,t) 
p(c(a,t),t) p(a,t = 0) da? 


Magnetic stress of the first order is: 

<5 (1) F m = |[<5 (1) B x rotB (d) ] + [B (d) x rot5 (1) B]| 

- ^-^[B (d) x rotB (d) ]. (22) 

47T/3 2 

The second order equations in a are: 

d t 8 (2) v + v k \7 k 8 (2) vi + v k X7 k 6 (2) Vi + 8 m v k V k 8 m vi 

= -Vi5 (2) U + + 5 {2) Fr, (23) 

d t S (2) p + v i X7 i 5 (2) p + S7i(p5 (2) v i ) = 0, (24) 

A<5 (2) <F = 4 n GS (2) p, (25) 


For the sake of simplicity we assume c l (a,t ) = a 1 + v l t, 
thus, stationary condition implies p(c(a,t),t) = p(a,t), i.e. 
the density p is constant along trajectory c(a,t). 

As the first step let us consider the drift solution. 
We define it as a perturbation of unperturbed station¬ 
ary configuration which is an exact solution of MHD 
equations. Let drift solution be given by the trajec¬ 
tory c^ d \t) = c(a,t) + ^ d \a,t ) and Eulerian variables 
(p^ d \r,t), f/ d ^(r, f), i/ d ^(r,f), B (d ^(r, £)). However, let us 
assume that Eulerian perturbations of <f>, p, p, v with respect 
to unperturbed stationary configuration, 8^ d 'f = f ^ — /, 
are second order in a: S^f = 0(a 2 ) (at least for a time 
of consideration) and their space and time derivatives have 
same order in a (i.e. they are not too large: < 

8^ f/R-, dt8^ f < 12<5^/).Here / is <F, p, p, v and a -C 1. 
We do not assume Eulerian perturbations of magnetic field 
to be small and thus do not expand them in series on a. 
Instead, we consider evolution of total magnetic field, which 
is given by the equation: 


p (d) (r, t) 


B i (a,t = 0) d(c (d) ) j (a,f) 
p(a, t = 0) da i 
B\a,t = 0) [,■ d(^y(a,t) 

p(a,t = 0) 1 da * 


■ ( 18 ) 


Let us now perturb drift solution by rapidly oscillating 
perturbation of amplitude a. To describe its evolution, we 
apply perturbed (with respect to the drift solution) MHD 
equations of the first and the second order in Eulerian form. 
In these equations one can substitute p^ d \ 4?^, and 
with correspondent values of unperturbed stationary 
configuration because they differ only in the second order 
in a. Thus, linear perturbation equations can be written in 
form: 


d t 8 w Vi + v k V k 8 (1) Vi + v k V k 8 w Vi 

d t S (1) p + uVi5 (1) p + V ; (p<5 (1 V) 
A<5 (1, <F 


-Vi8 (1) U + 8 (1) Fy, 

(19) 

0 , ( 20 ) 

4 nG8 (1) p. (21) 


10 Applicability of ideal MHD f or NS r-modes is justified by h igh 
conductivity of NS depths (see I Yakovlev &T Shalvbkovlll99lL for 
example). 

11 Strictly speaking, we require NS to be static in a rotating 
frame (corotating with star). We do not require magnetic field to 
be axially symmetric with the rotation a xis. We n eglect de forma- 
tions of a star due to magnetic stresses. lAsai et ahl (j 20151 ) show 
these deformations to be rather unimportant for r-modes. 


where second order magnetic stress is: 


5 2 F m = 


A- |[5 (2) B x rotB (d) ] + [B (d) x rot5 (2) B]j 


+ — \S (1) B x rot<5 (1) B] 

47rp 


<5 (1) p 
47rp 2 
8^ p 


j[<5 (1) B x rotB (d) ] + [B (d) x rot<5 (1) B] j 
(<5 (1) p) 


47rp 2 47rp 3 


[B (d) x rot B 


(d), 


(26) 


Here SB = 8 (1) B + S (2) B + ... = B (o) - B (d) is Eulerian 
variation of magnetic field on course of oscillations and B^ 
is total magnetic field (in accurate solution of MHD equa¬ 
tions). However, in contrast to J3 (d) , we are interested only 
to the first and second order perturbations of the magnetic 
field (with respect of drift solution), assuming Eulerian vari¬ 
ation SB ^ to be small. Let us present a trajectory of fluid 
elements for this solution in form 


c^(a,i) = c^ d '(a,i) + £ < ° , (c ( ' d ' ) (a,t),t). (27) 


This form allows us to use the same relations between Eule¬ 
rian variation of the velocity S^v and S^v and as we 
do in absence of drift motions [Eqs. <j9j) and (110M l 12 l 

Magnetic field and trajectory are coupled by the in¬ 
duction equation, which l eads to vanish ing La grangian 
perturbation of B/p (see IGlampedakis fe Aiiderssonll2007l . 
for example). Thus, Lagrangian variation A B l = 
—_B l Vk£^ fc and Eulerian variation of magnetic held can 


be derived 

following general formalism for vector 

helds 

(Friedman 

& Schuta 1978al:ISall2004h: 

8 w B (o)i 

— 

—B (d) i V 3 ^ 1)j -£ (1)j Vj B (d)i 



+ 


(28) 

5 (2) B (o) i 

= 

-B (d,i Vj £ (2)j ' + B (d)i Vj( (2)i 



- 

C (2)j Vj B (d)i - -£ (1)k £ (1)j \7 j 'V k B (d) 



+ 

V fe £ (1) i A (1) B (o) k - £ (1) k V k 8 (1) B\ 

(29) 


Equations (□MUI and mm differ from linear equa¬ 
tions in nonmagnetized case [mu and (00] only by mag¬ 
netic stresses 5 (1 ).F m and 8 (2 ^F m in Eqs. m and 123), 


12 Even the last equalities in Eqs. and EB can be applied 
for unperturbed stationary configuration corresponding to solid 
state rotation as far as drift motion perturb stationary velocity 
only in the second order in or. ~ a 2 v. 
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respectively. Let us analyse effects of magnetic field onto 
nonmagnetized ( B = 0) solutions by comparison of mag¬ 
netic and other stresses in these equations. Nonmagnetic 
(hydrodynamical) stresses can be estimated as « 

dtS^v = 0(a)Ruj 2 for the first order Eq. 11191 ) and as 
$( 2 )F h ss dtS^v = (D(a 2 )Rw 2 for the second order Eq. 
( 1231 ) . Before estimating magnetic stresses let us note, that for 
B = 0 solutions with nonvanishing drift, displacement (£°)^ 
contains terms oc a 2 cot, which are finite [i.e. (£°)^ = 0(1)] 
at tw > 1/a 2 . For such solution Eulerian variation S^B ~ 
Ba 2 Lot becomes large (5^B > B > 5^B). Presence of 
large Eulerian variation makes it doubtful, that the theory 
based on the expanding of Eulerian variations in series in a 
is appropriate here! 43 ! 

Thus, let us concentrate on the non-drifting solution 
at B = 0 (see Sec. EH) and demonstrate that it indeed can 
be applied to describe perturbations of the drift solutions 
(i.e. to describe oscillations at nonstationary background 
and B ^ 0). For this solution Lagrangian displacement 
stays always small = O(a) and even 0(a 2 ) terms do 
not increase with time, thus S^B = 0( a) ss aB^ and 
S {2) B = 0(a 2 ) ~ a 2 B ( ' A K Now we can estimate magnetic 
stresses of the first and the second order [given by Eqs. 
(HI and (25)] as 5 m « a(B w ) 2 /(4npR) and <5 (2) F m « 
q 2 (B^) 2 /(47t pR) respectively. Thus, magnetic to nonmag¬ 
netic stresses ratio in Eqs. (fl^l) and (l23l) is (5 (1) F’ m /5 (1) F' h « 
s (2) F ™/ 5 ( 2 ) F h fa (B< d ) ) 2 /(4irpR 2 uj 2 ) = v 2 a /(R 2 oj 2 ), where 
va = \J (B( d l) 2 /(47rp) being Alfven velocity. As a result, 
we conclude that magnetic stresses in both of these equa¬ 
tions are negligible, if B ^ <C B CIlt ~ Ruj^fFiTp. For a r- 
mode in NS and fiducial parameters (p = 4 x 10 14 gem -3 , 
w = 4/3 n = 8n/3u = 5000(iV600Hz) s -1 , R = 10 6 cm) 
we obtain B CI it ~ 2 x 10 17 G0 Thus, B = 0 non-drifting 
solution can be used as the first approximation to describe 
non-drifting oscillation mode for magnetized NS [first ap¬ 
proximation to the solution of Eqs. Ha and (|23ll25|) ]r^l 

Let us analyse corrections to this solution, which can be 
associated with magnetic field. In the first order these effect s 
were studied accura t ely bylLee (2005); Abbas si et al.1 (|2012l ); 
IChirenti &; Skakalal (|2013l )T TAsai et al.1 (|2015h . For simplic- 
ity, here we apply analytical estimates, which is enough for 
the purpose of the study. Magnetic stresses, which oscil¬ 
late along with oscillations, can modify oscillation frequency 
slightly by providing additional stiffness to oscillations. In 
analogy with effect of small magnetic field on sound waves 
in plasma, we anticipate increasing of the oscillation fre- 


13 It is worth to note, that large Lagrangian displacement is 
not a problem in nonmagnetized case, because Eulerian variations 
of all quantities are small at any moment of time and thus can 
be expanded in series in a. 

14 This estimate slightly differs from one obtained by 
iRezzolla et al.l d2000l) based on another physical arguments. How¬ 
ever, factor of few units difference is not a problem for order-of- 
magnitude estimates. 

15 This conclusion cannot be applied for non-magnetized solu¬ 
tions with nonvanishing drift (mass transfer velocity ^ 0, 

see Sec. n. Indeed, stress governing evolution of is zero in 

nonmagnetized case, but for magnetized star it should be com¬ 
pared with increasing [because ?4 ml ) jz o] magnetic stress F m . 
Thus, magnetic stress can not be ignored and should modify such 
solutions. 


quency to the value 5oj u> (v A /v s ) 2 , where v s is a phase 
velocity in absence of magnetic field. For r-modes, taken 
i jR as a phase velocity, we obtain Sui/w ~ (B^ /B C rit) 2 ■ 
This estimate is in a reasonab le agree ment with accurate 
calcul ation bv lAsai et al.l (l2015h and by IChirenti fe Skakalal 
(l2013ll (for rotation, which is not too slowll 16 ! The terms in 
S^F m with nonvanishing oscillation average (for example, 
oc B x rot S^ 1] B) cannot lead to significant displacement 
across the magnetic field, thus significantly affecting mag¬ 
netic field, for a simple reason: additional displacement of 
second order St; B = 0(a 2 ) across magnetic field would pro¬ 
duce compensating stress due to deformation of the mag¬ 
netic force lines. As a result, we suppose that such terms 
do not lead to significant modification of magnetic field 
and global drift motion. Concluding, for a given constant 
= const, none of magnetic terms leads to significant 
modification of nonmagnetized solution, which would affect 
evolution of B. In general case B^ is not a constant, but 
varies on a timescale ~ R/v^ ~ l/(a 2 uj) 1 /uj. 
Hence, adiabatic invariant of the first order oscillation equa¬ 
tion E/uj should be conserved during B^ evolution. As far 
as £ (d) < Bcrit, corrections to the oscillation frequency are 
small and oscillation energy is also conserved in leading or¬ 
der. 

Let us now demonstrate that a general second order so¬ 
lution can be presented in a form of superposition of drift so¬ 
lution and non-drifting r-mode. To do that, we analyse evo¬ 
lution of the oscillation mode with nonvanishing mass trans¬ 
fer velocity at the initial moment of time. It is espec ially im¬ 
portant to analyze this possibility as far as lFriedman et aD 
(l2015ll clearly demonstrated that r-modes excited by gravi¬ 
tational radiation in uniformaly rotating nonmagnetized NS 
correspond to well defined nonvanishing mass transfer ve¬ 
locity. Let us assume that at this moment Eulerian pertur¬ 
bations of the first order in a were given by an oscillation 
mode for initial magnetic field £T°Z In this case an arbitrary 
perturbations of the second order can be presented as a sum 
of second order perturbation, corresponding to the second 
order non-drifting solution and corrections, associated with 
nonvanishing mass transfer velocity. Let us use these correc¬ 
tions, which are second order in a, to define initial state and 
determine drifting solution c^ d ^(a,t) by solving exact MHD 
equations. As long as initial velocity perturbations are small 
[<5 < - d )v = u (mt) = 0(a 2 )], we expect that temporal evolution 
of this velocity field takes place on Alfven timescale, pre¬ 
serving Eulerian variation of all quantities (with possible 
exception of the magnetic field) to be of order of a 2 , i.e. tra¬ 
jectory c^(a,t) satisfies requirements, imposed to the drift 
motion solution in the beginning of section.Thus, we present 


16 For very slo w rotation (u < 0.07^/GM/R 3 ) 

IChirenti fc Skakalal j2013l ) suggests Suj/uj oc -B 4 . If this holds 
true, we can significantly overestimate the effe ct of the m a gnetic 
field on the oscillation frequency. However, lAsai et al.l (120151) 
suggest scaling 5uj/uj oc (B/B ct it) 2 (in agreement with our esti¬ 
mates) even for very slow rotation. The difference between the 
results of these papers can be associated with different magnetic 
field configuration and models of NSs applied by authors or, more 
probably, by too strong m agnetic fields B > 10 15 G discussed by 
Chirent i &; Skakala d2013l ). For such a large magnetic field the 
lowest order expansion 8uj/uj oc B 2 may become invalid at low 
rotation frequencies. 
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initial condition as a composition of drift solution and non¬ 
drifting mode. It allows us to apply this solution to describe 
evolution of initial perturbations. Thus secular evolution of 
the magnetic field follows the drift solution, i.e. it should be 
the same as in absence of non-drifting r-mode. As a result, 
we can estimate an upper limit for the enhancement of the 
magnetic field energy f» AB 2 R 3 /6 as an energy of initial 
drift motion « 2nR 3 p(5 (d) v) 2 /3. Here A B 2 = B 2 — (B m ) 2 . 
For initial drift velocity equal to mass transfer velocity 
= v (mt) « a 2 Ruj and fiducial NS parameters writ¬ 
ten above we come to %/A B 2 « 10 8 (a/10 -4 ) 2 G. Thus, 
if initial field B 10 8 (a/10 ~ 4 ) 2 G it will not be en¬ 
hanced by the drift significantlvl 17 ! If, on the contrary, ini¬ 
tial field B <C 10 8 (a/10 -4 ) 2 G, it can be enhanced up to 
10 8 (a/10 -4 ) 2 G. 

This section’s conclusion can be formulated in almost 
the same form as for nonmagnetized NS: a general second 
order r-mode solution can be presented as a superposition 
of oscillating solution with vanishing drift and a solution, 
which describes differential rotation in nonoscillating NS. 
The main effect of magnetic field is that it modifies differ¬ 
ential rotation in nonoscillating NSs. In case of magnetized 
star differential rotation cannot be stationary, but should 
evolve on the Alfven timescale. 


4 DISCUSSION AND CONCLUSION 

We reanalyze second order effects in r-mode oscillations. For 
nonmagnetized NS we refute widely believed statement that 
(second order) drift (secular motion of fluid elements) is an 
essential feature of r-modes. We present the second order 
r-mode solution with vanishing drift (non-drifting solution; 
see Sec. H. Strictly speaking, this solution is applicable in 
a case of slow rotation (12 <C 92k) and a (H/Hk) 2 - How¬ 
ever, we suppose that these restrictions are not crucial for 
the main results of paper (see footnote [ 8 ]), but we leave a 
strict proof of this conjecture beyond the scope of the paper. 
In Sec. [3] we demonstrate that this non-drifting solution is 
valid for magnetized NS. On the contrary, B = 0 r-mode 
solutions with nonvanishing drift cannot be applied directly 
to magnetized NS because magnetic field leads to evolution 
of the mass transfer velocity, which describes drift, on the 
Alfven timescale. However this evolution is decoupled from 
oscillations. This means that the general r-mode solution can 
be presented as superposition of two solutions: (a) the solu¬ 
tion characterizing evolution of mass transfer velocity (drift 
motion), which is the same as evolution of differential rota¬ 
tion in nonoscillating NS (initial conditions correspond to 
the mass transfer velocity at t = 0 ); and (b) non-drifting r- 
mode solution. Secular evolution of magnetic field is entirely 
determined by drift motion, thus its initial energy gives an 
upper limit for the enhancement of magnetic held energy. 
The total energy of the r-mode cannot be converted into 
magnetic energy because it exceeds drift energy for a fac¬ 
tor of 1/a 2 . Thus we conclude that magnetic held energy 
enhancement does not prove to be relevant for oscillation 
energy decrease and therefore cannot prevent instability. 

17 This estimate was suggested by M.E. Gusakov during our dis¬ 
cussions on the e arly stage of this wor k. I t is also in agreement 
with estimates bv lAlfven fe Felthammarl (Il963l l (see footnote n. 


Still, in our consideration we neglect gravit ational 
radiation-reaction force. As shown bv lFriedman et ah! (l2015h 
(published in arXiv, while the present manuscript has been 
under revision), this force excites r-mode with well dehned 
nonvanishing mass transfer velocity (differential rotation) in 
nonmagnetized NSs. Strictly speaking, in the present paper 
we demonstrate that this mode does not dump regardless 
of magnetic held en hancement (the opposite to what fol¬ 
lows from results by iRezzolla et al.l (l2QQ0l . l2001al lbll , in case 
when gravitational radiation force stop acting. In reality, 
gravitational radiation force acts permanently and accurate 
description of second order r-modes under combined action 
of gravitational radiation and magnetic held is particularly 
interesting task. We expect that two limiting cases are possi¬ 
ble, depending on the ratio of Alfven timescale 1 a to r-mode 
growing time r. If t\ <S r, then magnetic held has enough 
time to affect mass transfer rate and the theory described 
here remains applicable. In the opposite case, namely when 
tA t, r-modes will saturate rapidly, since they are almost 
unaffected by the magnetic held. 

Saturated r-modes require special consideration. It is 
typically assumed that saturation of r-mode instability is 
associated with the lowest parametric instability thresh¬ 
old in nonlin eary coup led trip l ets of oscillation modes 
[Schenk, Arras. Flanagan, Teukolskv fe Wassermanl l2002i; 
Arras. Flanagan . Morsink. Schenk. Teu kolskv fe W asscrman 


2003; 


Brinh J Teukolsk^ i 2^Wsserman 


[2004a b 


Bondarescu. Teukolskv fe Wassermanl [20071 . I200S 

Bondarescu fe Wassermanl 2013ll . Within this model 


saturated r-mode is always accompanied by excited daugh¬ 
ter modes (at least two). If amplitudes of these latter modes 
are comparable with r-mode amplitude, the hrst order ve¬ 
locity held should be presented as a sum of three oscillation 
modes, and each of them affects second order equations. 
Then the total mass transfer velocity, in principle, can differ 
from Eq. m, so that the very existence of nondrifting 
solution might be challenged as well as applicability of the 
theory discussed here. However, if saturation amplitude 
is not too large (a <C 1 ), we expect that such scenario 
does not occur and we can obtain nondrifting solution even 
for arbitrary number of eigenmodes excited in hrst order. 
However, a more detailed analysis of this question is left 
beyond the scope of this paper. 

Strictly speaking, the present results are accurate only 
at second order in mode amplitude. However, the drift 
of the fluid elements in the next orders (the third and 
higher) cannot be excluded. Even if mass transfer veloc¬ 
ity for such drift is very small v 3 ™^ « a 3 QR « 4 x 
10 -3 (iz/600Hz) (a/10 -4 ) 3 cms -1 , it can lead to large La- 
grangian displacements on a timescale of years (in nonmag¬ 
netic case). The second order drift and hrst order oscilla¬ 
tions are coupled in the third order. As a result the the¬ 
ory in third order can be much more complicated than the 
one in second order described here. Furthermore, for very 
large amplitudes a > 3, which can be studied only numer¬ 
ically, r-modes decay nonlinearly and strong differen tial ro¬ 
tation develops jGressman et al. 2002; Kast aunll201ll l. Such 
a strong differential rotation (angular vel ocities in the r ange 
0.5 ... 1.2 of the initial one according to lKastaunll201l} ) can 
strongly affect magnetic held (formally, it does not contra¬ 
dict with our estimates, which predict enhancement of the 
magnetic held up to ~ 10 17 G for a ~ 3). However, r- 
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mode s aturation amplitudes calculated bvlBondarescu et al.l 
(120091) ; iBondarescu fc Wassermanl (12013 1 are small a -C 1. 
Thus it seems reasonable to restrict consideration to small 
amplitudes. In this case, the total mass transfer velocity 
has small value of the second order in a , even if we in¬ 
clude high order terms. As a result, a general argument by 
lAlfven fc Felthammai] (119631 1 (see footnote [2|, that too slow 
differential rotation cannot generate magnetic field, is still 
valid within high order perturbation theory. Consequently, 
high order corrections can hardly result in strong magnetic 
field generation. However, we leave a more accurate analysis 
of high order effects outside the paper. 

To conclude we note, that results discussed here can be 
easily generalized for any oscillating mode in any magne¬ 
tized medium, if secular drift is a perturbation, allowed by 
the linear perturbation theory. This property seems to be 
rather general, at least it holds true for r-m odes and ocean 
waves IILonguet-HigginsI 1 19531 : iMoorel 1 197dl and we plan to 
study whether it is indeed a general property of second order 
oscillations or not. 
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